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Abstract 

The drawdown process Y = X — X of an one-dimensional regular 
diffusion process X is given by X reflected at its running maximum X. 
P*4 , The drawdown process Y — X — X_ is given by X reflected at its running 

(-H ■ minimum X_. In this paper we calculate the probability of the event 

{a a < (T|, A e,} for a, b > 0, where a a and at respectively denote the first 
passage times of Y and Y to a and 6, and e q is an exponential random 
variable which is independent of X. We then explicitly express the law 
of the occupation times of X below y until the first exit. Using this 
result we compute the laws of occupation times of X, its drawdown Y 
and drawup Y, at a a or e q . These results are applied to address problems 
{•f. , in risk analysis and for option pricing of the drawdown process. Finally, 

^^ ■ we present explicit formulas of the main results for the cases of Brownian 

f^"i | motion with drift and three-dimensional Bessel process, where we prove 

00 ■ an identity in law between an occupation time and the first passage time 

o 
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1 Introduction 



5_H . 

d . In a filtered probability space ($7, F, J 7 , P) with filtration F = {-?*}, we consider 

a one-dimensional regular time-homogenous diffusion X on / := (I, oo) with 
natural (or entrance) boundaries. Its evolution is governed by the stochastic 
differential equation 

dX t - n(Xt)dt + a(Xt)dBt, X = xel, (I) 

with an infinitesimal generator 

1 i, s ° 2 , s d 
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Here B is a standard F-Brownian motion, and (//(•), c(-)) is a pair of real- 
valued continuous functions. We introduce the running maximum and minimum 
processes of X by 

X t := sup X s and X t := inf X s , t > 0. 

0<s<t 0<s<t 

The drawdown process of X, and its dual, the drawup process, are then defined 
as Y := X — X and Y := X — X_, respectively. The first passage times of Y and 
Y above a positive threshold are respectively called the drawdown and drawup 
times. The occupation time of a stochastic process is the amount of time the 
stochastic process stays within a certain range. 

Laplace transforms of stopping times and occupation times for diffusion pro- 
cesses are well studied due to analytical tractability and their applications in 
risk theory, mathematical finance, and engineering. Some classical results for 
general diffusion processes using the Feynman-Kac representation and excursion 
theory can be found in [231 Hi]- Recent advances on this topic include jumps 
in the underlying model (see for example, [3J [T31 HI] ) , or incorporate memory 
into the model [7J. Among all the stopping times studied, one that is especially 
interesting is the drawdown time, which finds applications in financial risk man- 
agement (see for example, [U [HJ [T7J [27J 1323 ) , the theory of optimal stopping (see 
for example, [18, 22, 30 ) and the problem of the quickest change-point detection 
(see for example, jTOJ [T21 [23 [3TJ [35]). The Laplace transform of the drawdown 
time for general diffusion processes was first derived in [15]. Recently, [29] de- 
rived the probability that the minimum of the drawdown and drawup times 
proceeds an exponential random variable for a standard Brownian motion. In 
[37], the authors derived the joint Laplace transform of the drawdown time and 
the so-called speed of market crash for a general diffusion process, using pro- 
gressive enlargement of nitrations. Laplace transforms involving drawdowns for 
general spectrally negative Levy processes are studied, using excursion theory 
in [19] , which generalizes a result for diffusion- type processes in [TTJ EHJ [35l [36] . 

In this paper, we obtain a class of results regarding the sequential order of 
the drawdown and drawup times in a finite time horizon, and laws of occupa- 
tion times of X, its drawdown Y and drawup Y, for a general one-dimensional 
regular diffusion X. More specifically, we first derive the probability that the 
drawdown time precedes the drawup time in an exponential time horizon. We 
then compute the Laplace transform of the occupation time of X below y until 
the first exit time. Using this result, we proceed to compute the Laplace trans- 
form of occupation times of X below y, of Y above y, and of Y below y, until 
the drawdown time or an exponential time independent of X. 

Rather than relying on standard techniques in excursion theory, we use a 
perturbation approach to obtain the Laplace transform results. To the author's 
knowledge, the first use of this approach dates back to [15], where the condi- 
tional Laplace transform of the drawdown time is obtained by an approximation 
argument based on the Laplace transforms for the first exit times. Thanks to 
the classical theory of second order ordinary differential equations and Ito's dif- 
fusion, the Laplace transforms of the first exit times can be expressed in terms 



of eigenfunctions of the underlying diffusion process. Recent applications of this 
approach can be seen in [13] for refracted Levy process and [16] for diffusion 
process. 

The results obtained can be applied in risk analysis and for option pricing of 
the drawdown processes. In particular, we consider a time-homogenous diffusion 
with reduced form default model. The probability of realizing a drawdown 
before a drawup before the default time can be computed when the hazard rate 
is a constant. Using the so called Omega model (see for example, [T] [5]), we 
can describe the hazard rate of the default in such a way that it depends on 
the asset process, its drawdown or drawup processes. Then the probability of 
default before a large drawdown can be computed. Moreover, our results can 
be used to price Parisian-like digital call options and a-quantile options of the 
drawdown process, a non-trivial extension of the option pricing problem for 
maximum drawdowns g] US] [35] [33] [3JJ] . 

As examples of our general result, we present explicit formulas for some of the 
main results in the cases of Brownian motions with drift and three-dimensional 
Bessel processes. Moreover, we prove through Laplace transform that, in these 
two models, the law of the occupation time of the drawdown process above a 
level is the same as that of the drawdown time of certain threshold. 

The paper is structured in the following way: standard identities regarding 
the first hitting time and the first exit time of an one-dimensional regular dif- 
fusion are reviewed in Section 2. In Section 3, we derive the probability that a 
drawdown precedes a drawup in a finite time-horizon. In Section 4, we use the 
perturbation method to derive the Laplace transform of the occupation time of 
X below a level until the first exit time. We then use this result to compute var- 
ious occupation time formulas regarding the drawdown and drawup processes. 
In Section 5, we discuss applications in risk analysis and for option pricing of 
the results obtained in Sections 3 and 4. In Section 6, we present explicit for- 
mulas for some of the main results in the cases of Brownian motion with drift 
and three-dimensional Bessel process. We conclude with some closing remarks 
in Section 7. The proofs of propositions and theorems omitted can be found in 
Section [S] 

2 Preliminaries 

Let X be the linear diffusion process defined in (jTJ) . The first hitting time and 
the first passage times to a level y G / by X are respectively given by 



T, 



± 



inf{t > : Y t = y}, 
= inf{t > : Y t = y}. 



v ' L ' z < 



Throughout the paper, we use P^(-) to denote the measure P(-|.X"o = %), and E^ 
the expectation under ¥ x . Moreover, for q > 0, we denote by e q an exponential 
random variablq^J with parameter q, which is independent of X. 



As a convention, we assume that P(eo = oo) = 1. 



It is well-known that, for any q > 0, the Sturm- Liouville equation (Cxf)(x) = 
qf(x) has a positive increasing solution </>+(•) (decreasing solution (j)g(-), resp.). 
In fact, for an arbitrary fixed y 6 /, we can choose 

(E x {e-^}, ifx<y ( 1 if 

<t>t(x)={ 1 , cf>-(x)=\®y{e-^y - y ,VxeI. 

U^e-^}' ^ [E4e- ?T «}, i£x>y 

(2) 
The functions </>^(-) can be used to define: 



W q (x,y) := dct 
with its derivatives 



g («) </> q (y). 



Vs.i/ez, (3) 



W,,i(a:,i/) := ^^(x,y), W,, a (a;,i/) := g-W qil (x,y). (4) 

Without loss of generality, we assume that <j)j{x) satisfies (see, for example, 
page 19 of 0) 

W q ,i(x,x) =s'(x), Vx€l, (5) 

where s(-) is a fixed increasing scale function of X. When q = 0, we extend the 
definition of W q using 

W (x,y):=s(x)-s(y). (6) 

And we assume that W q {x,y) and its partial derivatives are continuous in q at 
0. The functions W q and 4>~ have the following properties: 

Lemma 2.1 For any x,y, z 6 I, q > 

wifcrt—wifcx), s^ - 4g^.&»>. jst*rw— . 

We recall the following result regarding the first exit of X from page 603 of [TH] . 



Lemma 2.2 Suppose that x,y, z £ I and (x — y){z — x) > 0, for q > 0, we have 
E x {e~^ ;Ty<Tz] = ¥x{Ty<TzAeq) = ^|^j. (7) 



3 Drawdowns and drawups 

In this section we study the law of drawdowns and drawups of X, where X is 
the linear diffusion process defined in ([1]) . To this end, we introduce the running 
maximum and minimum processes of X by 

X t := sup X s and X_ t :— inf X Sl t > 0. 

o<s<t 0<s<t 

The the drawdown and drawup processes of X are defined respectively as 

Y t := X t ~ X t and Y t := X t - X t , t > 0, 

and the drawdown of a units and the drawup of b units are the following first 
passage times: 

a a := inf{t > 0, : Y t >a},a> 0, (8) 

a b := inf{t > : Y t > b}, b > 0. (9) 

In the remainder of this section, we will derive the probability that a a pro- 
ceeds &b in a finite time interval [0, t] for any t > 0. This is a nontrivial extension 
of the infinite time-horizon result in 28 , which enables us to study the draw- 
downs and drawups for diffusions with killing, as well as limiting conditional 
distributions of a a as the thresholds a, b tend to infinity in a property manner. 
The task is accomplished by computing the Laplace transform of the function 
Px-.a.b(t) '■— ^"x(ca < &b A t) . In particular, for any x,a such that x — a G /, we 
will calculate: 

/>oo 

P*(ct q <a b Ae q )=q e~ qt P x (a a <a b At) dt. 
Jo 

3.1 The case of a > b > 

Theorem 3.1 On the event {a a < <7& A e q }, we have cr b < a b and X ab 6 

(x — b, x). Moreover, for u 6 (x — b, x) we have 

^x(o- a <a b A e q , X aa e b - a + du) 

_ s '(u + b)w q (x,ui ex fr Wi(v,v + b)\ du _ 



W^(u + b,u) \Ju+b-a W q (v,v + b) 

On the event {a a < o b A e q }, we have a b < a b and X& b G (x,x + b). Moreover, 
for u £ (x, x + b) we have, 

V x (a a <<r b Ae q ,X aa ea-b+ du) 

_s\u-b)w q{ u, X )f_ r +a - b ^v- b ) dv ) du . (11) 



W*(u,u-b) "V L W q (v,v-b) 



Proof. First, it is easily seen that for t > 0, and u € (x — b, x), 

{a b e dt,X t e du,a b > t} = {t u E dt,X t 6 u + d&}, P x -a.S. 
It follows from integration by parts that, 

^x{<y b <ffiAe,; X ab E du) = / ge _9f P* (ct 6 < a b A i, X CTb e du) di 

Jo 

/>oo 

= / c- qt P x (a b e dt,X t e dii,o- 6 > t) 
Jo 

/■OO Q 

= ^ e^-P x (r-edt,r+ +fe >t)dn. 



From Lemmas 12.21 and 12.11 we obtain (|10[) for the case a = b. 

If a > 6, then any path in the event {er a < a b A e q } can be decomposed into 
two path fragments: {X t }o<t<a b and {X t } ab <t<a a - Conditioning on {X ab =u}, 
the second path fragment is a process starting at u, and decreasing to u + b — a 
before it incurs the drawup of b units. Formally, using Markov shifting operator 
(X t o 6(s) = X t+S ), we have 

(?a= <Jb+ T^-a+b ° ^b), Vm £ (x - 6, x) . 

Using strong Markov property and memory less of e 9 , we have 

F x (a a <a b A e q , X aa e b - a + du) 

= V x {a b <a b A e q ,X ab e du} ■ F u (r~_ a+b < a b A e q ). 

Now ([TUl) follows from Proposition 13.21 below. Eq. (fTTj) can be proved using a 
similar argument. ■ 

Proposition 3.2 For x,y (z I, we define m — max(i, !/) and n = miri(x,y). 
Then we have 

'-<*<**-> - •*(f 3 §g&w4 <12 » 

*.(«,. A „) = ^(-jT ^^ d.). (13, 



3.2 The case of b > a > 

To obtain the result for general b > a, we recall the following result for the 
limiting case b — oo from page 602 of [15] . 

Proposition 3.3 

,.(„.<„,_ r 'w expf- ri|i<-i^) d „) d „. (14) 

A Wq(u,u-a) V A Vv g (v,v-a) / 



A useful observation is that 

P x (va < o-h A e g ) = P x {a a < e q ) - ¥ x (a b < a a < e q ). (15) 

The second term on the right hand side of the above equation is computed in 
the following result. 

Theorem 3.4 On the event {a b < a a < e q }, X„ a G (x,x + a). Moreover, for 
u G (x, x + a) we have 

V x {&b < o-g < gg, X da £b- a + du) 
du 
= s'(u-a)W q (u,x) t _ f u+b - a W qA (v,v~a) ^ 

Wjj(u,u-a) V Ju Wq(v,v-a) 

s'(v) _( f v W qA (w,w-a) 



exp — / — ^— — ■ —dw dv. (16) 

lu+b- a W q {v,v-a) \ J u+b _ a W q (w,w-a) J 

Proof. Using strong Markov property and memory less property of e g , we have 

¥ x (a b < a a < e q , X„ a £ b - a + du) 

= P x {a b < e q ,X aa e b - a + du) •P M+b - a (c r a < e s ). 

The result follows from Theorem 13. II and Proposition [ 



Remark 3.1 By letting q — > + in M0\) . Hi]) , (ITSl) and \16\) , and using ([6]), 
we obtain Theorems 4-1 o,nd 1^.2 in '281. Second, by letting b — y oo in (|15|) 
and Iil6\) . one obtains the Laplace transform of a a . The results in Theorems 
\3.1\ and \3.J\ also enables us to compute E x {CTa|<T a < crt>} and the limit law of o~ a 
given o~ a < a b , as a, b —¥ oo in certain way. These distributional results can be 
applied to the problem of sequential detections and identification, which we leave 
as future work. 

4 Occupation time formulas 

In this section, we begin by computing the Laplace transforms of the occupation 
time below a level until the first exit time for a linear diffusion process X defined 
in (p}. Although relevant formulas exist for special diffusions, the results for 
general linear diffusions that we provide here are new. Using this result, we 
then proceed to study occupation times of X, its drawdown Y and its drawup Y 
until the drawdown time a a or until an exponential time which is independent of 
X. These new results give several interesting fluctuation identities and provide 
means of measuring financial risk and pricing options as we shall in the next 
section. 



4.1 Occupation time below a level until the first exit time 

In this subsection we study the law of occupation time until the first exit time. 
In particular, for y S (a, b) C I, the occupation time below y before exiting 
(a, b) is denoted by 



A a v b := l{x t<y} dt. (17) 

Jo 

The law of Ay ,b is summarized in the following result, the proof of which can 
be found in Section |8l 

Proposition 4.1 For x, y € (a, b) C J, q > 0,]? > 0, we ftaue 



E x {e-^ , -^;r+<r-} 

'iy ?+p (x,a) s'(y) 



= < 



wr ^ (y>o) ^i(v,6) + iy,(6,v) ^ (y 'f 

W p (x,y) , W p (b,x) s'(y) 



W P (b,y) W p (b, y ) WpAy>b) + WpM W q+Ptl (y,a) 



W q+P (y,a) 



E x {e-« A Z' b ;T+>T-} 

W q (y,x) W q (x,a) s'(y) 



W q (y,a) W q (y,a) xxr ^ _, , s'(y) 



W q Ay,a)+ s{b) _ s{y) W q (y,a) 



(s(b)-s(x))s'(y) 



I (s(6) - s(y))W qtl (y, a) + s'(y)W q (y, a) ' 



ifxe (a, y] 

*/ * G (y, b) 

(18) 

i/ a; e (a, y] 

*/ a: G (y, b) 

(19) 



Remark 4.1 Lemma \2.2\ can be recovered from Proposition \4-l\ by the limit 
y —} o— in (|18p and (|19| . Moreover, Corollary 3.1 in \16^ can be easily obtained 
by letting x = y, b — > oo and a — > /+. 

Letting a — > /+ in (fTHJ) and using Lemma ETT1 we obtain the following result. 



Corollary 4.2 For x G (y, o) C /, p > 0, q > 0, we ftave 

E x {exp i-q Mx t < v } di - pr b + J } 

W p (z,y) W p (6,a:) s'(y) 



W p (b,y) W p (b,y) rir f jy rrrn Jt+p'W 



W pA (y,b) + W p (b,y)^ 



J q+p 



(y) 



4.2 Occupation time below a level until the drawdown 
time 



For 1-0,1/eI, the occupation time below y before the drawdown process Y 
hits a is denoted by 

B a y := ri {Xt<y} dt. (20) 

Jo 

While the occupation time of X below y until an exponential random variable 
independent of X is well-studied in literature [16] , the new quantity By defined 
as above relates the occupation time of X to its drawdown process Y, which can 
be used to characterize the drawdown risk of X. The law of B° is summarized 
in the following results, the proof of which can be found in Section [3] 

Theorem 4.3 For q > and x G I, 

e x {e-* B * } = exp f - / — : {x ; u :t , r 

V X 1+ s'(x) W q (x,u-a) 

x+a 1 s'(u)du 



(, , s(u)-s(x) WVi(:r,«-a) \ W q (X,U - a) 
y 1 "I V(x) W q (x,u-a) ) 

4.3 Occupation time of the drawdown process until the 
drawdown time 

For any y £ (0, a) and x — a € I, the occupation time of the drawdown process 
Y above y before Y hits a is denoted by 



Cy : = / hY t>y} dt. (21) 

Jo 

The occupation time C° measures the amount of time for the drawdown process 
Y to finish the "last trip" from y to a. It can be used as a measurement of 
performance for CUSUM-type stopping rule in change-point detection problems 
[2"5] . Because of its obvious financial interpretation, C" can also be used as a 
measure for drawdown risks. 

To obtain the law of C", we first condition on X aa , and then count separately 
the occupation time before and after the moment when the peak X aa is realized. 
To this end, we recall the following result from page 602 of [15j . 



Lemma 4.4 For m > 



x. 



.■(A' „ :;/,;) = <-.np( - / S ' {V ) -dv). (22) 

' s(u) — s(v — a) 



Theorem 4.5 For q>0, < y < a, 



E x {e 



-*c; 



a a < 00} 



s (m) 

VKg (m — y,m — a) 
-. , s(m) — s(m — y) Wq,x(rn—y,m— a) 

s' (m — y) W q (m — y.rn — a) 

s' (u) W q i(u—y f u—a) 

,_, — ppj^ppr 

J s(uj — s(u — y) Wq,i(u—y,u — a) 



1 i ■=> V " ) — &yui — y j vv q.iyu, — y,u — u, 
s' (u — y) W q (u—y,u — a) 



dm. 



Proof. Following the conditioning argument in [37] , we define 

9a ■= sup{i < a a : X t = X t }. (23) 

Given {X 9a = to}, the path fragments {^t}te[o,g a ] an d {Xt}t£[g a ,a a ] are two 
independent conditional processes. For any to > x, we let e = "'^ for a large 
integer N > 0. Using Lebesgue dominated convergence theorem, continuity and 
strong Markov property of X, we have that 



E^jexp f - q / l{y t>1/ } di j | X 9a = m} 
E K {expf-q/ l { y f>y} dt) |cr a >r+} 

iV-l 

lim TT E x+ie {e 



A x+U- a ,a: + (i + l) t 



J 



JV-1 



exp hm N \E x+ie {e- qA »+™-v 



T x + (i+l)t ^ T a;+ie- 



exp 



*'(«) 



s'(u) W q ,i (u—y,u — a) 
s' (u—y) W q (u — y,u—a) 

s(u) — s(u — a) I _|_ ^(^)-s(u-2/) W q ,i(n-ij,u-a) 

s' (u — y) W q (u — y,u — a) 



d«), 



where we used Proposition ^. ll in the last equality. Similarly, for the occupation 
time after the random time g a , we have that 



Ejcjexp - q / l{y t>y } di | X 9a = m} 



.{exp [ -q l{x 4 <m-2/} di ) | r m _ a < r+} 

l+e' 



lim ■ 

e'->0+ 



s{e 



"9^„ 



) 'm— a ^ 'm+c'J 



H T m-a < T m+e'l 



s(m) — sirn — a) 
W q (m — y,m— a) 



-. . s(m) — s(m— y) W Qs i (m— y,m— a) 
s' (m — y) VKq (m—y,m—a) 



The result now follows from Lemma [44 



10 



4.4 Occupation time of the drawup process until the draw- 
down time 

For any y £ [a, oo) and x — a £ I, the occupation time of the drawup process Y 
below y before the drawdown process Y hits a is denoted by 



D v--= / h%<v} dt - ( 24 ) 

J 

The occupation time D y can be considered as a counterpart of C" defined in 
(j2"Tj) . The law of D" is summarized in the following result. 

Theorem 4.6 For all q > and y > a, 

rx+a s'(u-a)W q (u,x) 



E x {e- qD y ; o- a < oo} = P x (a Q <a„Ae,)| / 

J X 



W^(u,u — a) 



V a Wg(v,v-a) ; 

Proof. Notice that D° = er a , P^-a.s. on the event {a a < a y }. On the other 
hand, on the event {a y < a a }, we have X_a — X& y — y < X& — a < X Ga — a, 
P x -a.s. Thus, 

D a y =a y + B a Xay o9(a y ). 

Using strong Markov property we have, 

E x {e- qD v;a a <oo} 
= E x {e~^;a a < a v } +E x {e-^l {ay<aa} E x , v {e' qB K}} 

i-x+a 

= Px(o-a < o y A e 9 ) + / P x (v y <o- a A e q , X dy ey - a + du) 

J X 

xE„_ a+9 {c-« B "-»+»}. 
The result now follows from Theorem 13. II ■ 

4.5 Occupation time of the drawdown process at an inde- 
pendent exponential time 

For any y G (0, oo) and x — y € I, the occupation time of the drawdown process 
Y above y before an independent exponential time e q is denoted by 

E* := r q i {Yt>y} dt. (25) 

Jo 

The occupation time E^ can be identified with the Laplace transform of the 
occupation time in a finite time horizon T > 0: 



Mv^yydt, 



11 



which is closely related to the maximum drawdown sup te r 0]T i Y t and other quan- 
tiles of Y in the finite time-horizon T. The law of E% is summarized in the 
following result, the proof of which can be found in Section [5] 

Theorem 4.7 For all q,p > and y > 0, 

^x{c p «} = l-cxp - / ^±P — — du 

V h W qil (u ~y,u) + W q { Ul u-y) fg^g 

XP ' ' ^> + '(u-y) 



Q + P 



W qA (u -y,u) + W q (u, u-y)— 
-E-a / rm) « + '' ,(m ~ y) 

g+P l j < +p (m^) 

W q: i(m — y,m) + W q (m,, m — y)- 2±£ 



dm. (26) 



5 Application 

5.1 Probabilities regarding drawdowns and defaults 

A realization of a large drawdown is usually considered to be a sign of market 
recession. In this section, we use a reduced form model for default and compute 
the probabilities of drawdowns and default. In particular, we consider an asset 
process X, which is given by a time- homogeneous diffusion process with initial 
value x and lifetime £: 

dX t = fi(X t )ddt + a{X t )dW t , t < C (27) 

Here £ is an independent positive random variable which models the "default 
time" of asset process X. If we assume that the "default time" £ = e q for a 
q > 0, then the probability that there is a drawdown of a units before a drawup 
of b units by the default time is given by ¥ x (<7 a < a^ A e q ), which is readily 
available from Theorems 13.11 and 13.41 

Moreover, we can consider a more realistic model for the default time C,, to 
reflect the fact that realizations of drawdowns of the asset process are very likely 
to be followed by a default. In particular, we adopt the Omega model studied 
in [TJ |5J [T3] to model the hazard rate of £ at time t > as ql^ Yt > y } for some 
y € (0,a). Here a is a large number that characterize a critical level of a large 
drawdown. Then the probability of default before the drawdown of a units is 
given by 

P*(C < O = 1 - P,(C > <r a ) = 1 - K{e- gC y}. (28) 



The Laplace transform in (J28J) can be found in Theorem 14.51 

Hazard rate of similar form can be considered. For example, we can model 
the hazard rate of the default time £ as ql{x t <x}- Here the initial value a; is a 
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critical benchmark level which may trigger a default through default intensity. 
We can also model the hazard rate of the default time C as <fi-s Y < v \- This 
is the case in which the default tend to occur when there is not enough upside 
momentum for the asset process. In both cases, the probability of default before 
a drawdown of a units can be found using Theorems 14.31 and [ 



5.2 Option pricing for the drawdown process 

Options on maximum drawdown and drawdown processes have drawn lots of 
attentions in recent years (see [4j[5j[26, 32, 34, 39 ). In this section, we use an 
semi-analytic approach to price a large class of options on the drawdown process. 
In particular, we assume that the market is complete, P is the risk-neutral 
measure, and r > is the risk-free interest rate. We model the underlying 
procescl as the time-homogeneous diffusion X defined in (TTJ). Then a Parisian- 
like (see [6] for a definition of standard Parisian option) digital call on the 
drawdown process with barrier y > 0, maturity T > and strike K £ (0,T) is 
worth 



/ Mv t >y}dt>K), 
Jo 



P (x,y,K,T) = e~ rI F x ( l {Yt>y} 6t > K), (29) 



at its inception. Using double randomizations: T — e q and K = e p , then the 
randomized option price is given by 

P (x,y,q,p) = e~ rT - e- rT E x {e~ pE y}. (30) 

The Laplace transform in (|30p is readily available in Theorem 14.71 Hence, the 
price ((29)) can be computed via double Laplace inversion: 



P (x,y,K,T) = e' rT - e- rT -l p (- -iJ-E^e-^y} 

\P \Q 



(31) 



K 



where X p and I q are Laplace inversion operators. 

Corridor options such as a-quantile option are studied in [20j [21] . Below we 
consider pricing of a-quantile options on the drawdown process. To this end, 
for anae (0, 1], we define the a-quantile of the drawdown process during [0, T] 

by 

Y« := in% > : / l {Yt>v} dt < (1 - a)T.} (32) 

Jo 

In particular, Y^ = sup te r 0T i Y t is the maximum drawdown at time T. An 
option on the a-quantile with maturity T > and an absolute continuous, 
bounded payoff function /(■) such that /(0) = is worth 

A (x,f,T) = e- rT E x {f(Y£)}, (33) 



2 Notice that the underlying process is not necessarily an asset price process. It can be, for 
example, the logarithm of an asset price process. 
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at its inception. We notice that 

M/OtO} = E A l { Y°>u}f'(u)du}= f'(u)-P x (Y«>u)du 

/'(«) • P*( / l{Y t >u}dt > (1 - a)T) du. 
/o Jo 

It follows from ^ that 

A (x,f,T)= f'(u)-P (x,u,(l-a)T,T)du. (34) 

Jo 

Again, by double Laplace inversion and Theorem 14 . 71 we can compute the price 
in dMJ). 

6 Examples 

6.1 Brownian motion with drift 

In this section we derive a group of explicit formulas for a Brownian motion 
with drift. In particular, we consider a Brownian motion with drift /i ^ and 
diffusion coefficient a > 0: 

dX t = ji dt + adB t , I = (—00,00). 

Let us denote by 

Then the increasing and the decreasing eigenfunctions of X can be chosen as 
(see for example, [5]) 



It follows that 



sM^(l-e- 2fe ), 

u/ c \ o -*(*+v) sinh[7(o; - y)] Wg.ifoy) 

W,(x,y) = 2e °^ + ^ , — 7- =jcosh[j(x-y )} - S. 

1 Wq{x,y) 



From Theorems 13.11 and 13.41 we have that 
Corollary 6.1 

PoK <o-& Ae,) 

Sb 



V 2 7 / e Sb ( 1C Qth[nfb]+S) _ 7 \ e (q-6)(5- 7 co S hhfe]) > j, Q . 

2g V sinh[ 7 6] sinh 2 [ 7 6]y 



I cr 2 7 / G^ a (7coth[7a]-i5) 7 \ c fa-b')([5+7 coshhal) 



2q \ sinh[7a] 

-fe~ 5a , b > a > 0. 



7 cosh [7a] — <5sinh[7a] 
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o(o-& < <J a Ae,) 

' ^1 f ^ a (7COth[ la }-S) _ 7 V(a-fe)(5 +7 c OS h[ 7 a]) i & > a > ; 

2q \ sinh[7a] sinh 2 [7a]/ 



<t 2 7 f e l5b (7coth[7b]+<5) _ 7 ^ p (6-a)(-o+ 7 cosh[7b]) 

«6 



1 _ g 7 g i,7<juiii[7oj-t-u,) 7 I „ 



2g ^ sinh[7fc] sinh 2 [76] 

7e 00 , a > & > 0. 



7 cosh [7 6] + (5sinh[7&] 
From Theorem 14.31 we have that 



Corollary 6.2 

w f - qBS 9 . , f f e*" 6( 1 coth[ 1 (a^u)}-6)du 

E {e q °;X CTa >a\ = exp - / . r , -rT7 ^TT^ — Tul - 

\ J sinh|du| 7 cot h [7 (a — u)\ + o coth|du 

E {e <?±i o ; X CTa <a} = / '■ 1—* J- 

Jo sinh[7(a — u)J (7Coth[7(a — u)\ + 5coth[Ju]) 



For occupation time of the drawdown process, from Theorem 14.51 we have 
that 

Corollary 6.3 For any y € (0,a), the occupation time C c y L has the same law as 
the drawdown time o~ a -y 

Proof. Straightforward calculation yields that Using Theorem 14.51 we obtain 
that 

s'(m) 



Eo{ e -9°?} = f" ii.,,,„- w .,„-„, 

Jo 



, s{m)—s(m — y) W qt i(m-y,m-a) 

s' (m—y) Wq(m—y,m—a) 

/ r m s'(m) ty g ,i(ti-y,«-a) 

X exp | - / ^'(u-a) W^-y.u-q) \ dw 

V /n 1 J- s(n)-s(tt-y) W^i (u-j/.M-q) / 

V U "*■ ^ s'(u-y) W q (u-yM~a) 7 

7 cosh[7(a — y)] — 6 sinh[7(a — y)] 



where the last equality follows from Proposition 13 .31 or Corollary 16. II as b — Y oo. 
It follows that the occupation time C° has the same distribution as (J a -y under 



6.2 Three-dimensional Bessel process (BES(3)) 

In this section we study the case of three-dimensional Bessel process. In partic- 
ular, we consider 

dX t = — dt + dBt, 7 = (0,oo). 

Let us denote by 

v:= V^q. 
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Then the increasing and the decreasing eigenfunctions of X can be chosen as 
(see for example, [2]) 



It follows that 



, + , s s\nh[vx] 



-, W^(«,j/) = smh i/(a; - y) , — f- - = — 

x vxy W q {x,y) x 



s(x) = --, W q (x,y) = 7:zr smh[u(x - y)], Jr r„\.\ = - - + ^ cotb.[i/(a; - y)] . 



Using Theorem 14.51 we have that 

Corollary 6.4 For a; > a > y > 0, the law of the occupation time C" is the 
same as the drawdown time <J a -y 



Proof. Straightforward calculation yields that Using Theorem 14.51 we obtain 
that 



s' (in) 

j'^JKJ i '. ^- 

^-qC^\ _ j W q {m-y,m-a) 

s('tn) — s(m — y) W qj i{m — y,m — a) 
s' (m—y) W q (m—y,m—a) 

s (u) VK Q) i (u—y,u—a) 

{ s'(u-y) W q {u~y,u-a) 

X GXp I 



/•OO 

Jx 1 



cosh[i/(a — y)] 



-, , s(u) — s(u—y) W q ,i(u—y,u-a) 
s'(u — y) W q {u — y,u—a) 

x-(a-y) , tanh[Ka- y )]' 1= ^ {c _ , , ^ 



where the last equality is obtained by substitutions a — > a — y and y — > 0+ in 
(|53"|) . It follows that the law of C° is the same as that of a a - y . ■ 

Remark 6.1 The results of Theorems \6.3\ and \6.J\ show an nontrivial fact: if 
X is a drifted Brownian motion or a three-dimensional Bessel process, then for 
a fixed y > 0, the law of a v is the same as C% +a for any a > 0. That is, the 
amount of time the drawdown process Y spends in [a, a + y] until the drawdown 
time a a +y is the same as the drawdown time o~ v . 

7 Conclusion and future work 

The perturbation approach is a powerful tool to reduce Laplace transforms of 
hitting times involving the running maximum/minimum, or Laplace transforms 
of occupation times, to that of the first exit time. In this paper, we apply this 
approach to derive the probability that a drawdown proceeds a drawup in an 
exponential horizon under a regular diffusion process X. We then proceed to 
compute the law of the occupation time of X below a level at the first exit time. 
Based on this result, we explicitly express the laws of the occupation time of X, 
its drawdown process, and its drawup process, until the drawdown time or an 
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independent exponential time. These results are applied to two problems in risk 
analysis and option pricing of the drawdown process. As examples of the general 
results, we present explicit formulas for some of the main results in the cases of 
Brownian motions with drift and three-dimensional Bessel processes, where we 
also prove a nontrivial identity in law between the occupation of drawdown and 
the drawdown time. An interesting extension of this work is to allow jumps in 
the underlying model, such as considering a spectrally negative Levy process or 
a jump diffusion process, and we leave it as future work. 

8 Proofs 

Proof of Proposition [3721 First we notice that 

P m (r~ <a b Ae q )= E m {e" 9T " ;t,7 < a b }. 

To compute the above expectation on the right hand side, we follow the idea 
of [15j . and partition the interval [n, to] into N equalength subintervals with 
length e = m ^ 21 - In particular, 

E m {e-« T ";r,; < a b } 
= Km E m {e^ E -» 1(T ™-P+i)."''"™ ) ;r ,.,,, < r+ , , ., VO < j < N - 1} 

N-l 

= jfe, II E™- 4e {e- A —(* +1 ). ; r"_ (i+1)6 < r+_ ie+ J. (36) 

i=0 

Here we used the Lebesgue dominated convergence theorem in the first equality, 
the strong Markov property of X and continuity of the running minimum process 
X in the second equality. 

To compute the limit in (|36j) . we use Lemma [2721 to obtain that 



N . . 

4=0 



N-l 

X 
e ^ ' m—ie+b i 



Urn J] E m _ ie {e ^w* ;r m _ (i+1 < r,i 



i. TT W„(m — ie,m — ie + b) 
lim 



JV-i-oo J- 1 W q (m — (i + l)e, m — ie + b) 

, N-l 

exp lim > 



W q {m — ie, to — ie + &) 



JV->oo 

i=0 



W q {m — (i+ l)e, to — ie + 6) 



ex P / -T7T7 JT" du 

1 /„ W g (u,u + 6) 



which completes the proof of (|T2j) . Eq. (1131) can be proved using a similar 
argument. ■ 
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Proof of Proposition ^, ll We follow the perturbation method in |T3] [TB]. 

For e > such that y + e < b, we approximate A y ' b by A y ' h e : 

OO 

A vi := z2( T y+t A Tq ^ ~ T y' n A r °> b )> 

n=l 

where t^' 1 := t~, and for n > 1, 

r+£ := inf {t > r"'" : X, > „ + e}, v" +1 := inf{i > r+f e : X t < y}. 
Using strong Markov property and the continuity of X, we have that 

E y {e-« A y*-ri;T+<T-} 
= Ey{e-^K + .. T ^ < r-} ■ E y+e {e-^:-^ + ;r+ < r"} 
= E,{e-^'< + .;r+. £ < r"} • (E. H+e { e -^ + ; r+ < r"} 

+E, +e {e-^ ;r - < T + }Ey{c - q A^ P r+. T + < T - } ^ 

from which we obtain that, 

E, +e {e-^ + ;r+ < r-}E,{c-^K + + .; r+ +e < r~} 
1 - E y+e {e-^ ; Ty < r+}E y {e-^<^ ; r+ +e < r a -} 

„. , - Ty g+P fa,a) 

WpU/ + £ 1 ^j W g+p fa+e,a) 

Tyjb w) i W p (6,;/+6) W, +P (y,a) 

PK "' ^ W p {b,y) W q+p (y+e,a) 

Using Lebesgue dominated convergence theorem and the continuity of X, we 
have 

E y {e-^-rf ; r+ < r~ } = ljm E y {e^ A ^'^ ; r+ < T ~ } 



It follows that, for a; € (a,y), using strong Markov property of X, we have 

E x {e-^' 6 -^V+<T-} 
= E x {c-^ + ^y;r+ < r~} ■ E y {e^ A y b '^; r+ < r"} 
W g+P (x,a) y(y) 

W q+p (y,a) WpAyib) + Wp ^ y) ^M- 
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For x G (y, 6), we similarly have 

E x {e-« Aa » b -ri;r+<T-} 
= E x {e~^ + ; r+ < t- } + E^e-^v" ; T - < T +) . E y {e"*^ 



' 'b 



Wpfoj/) , Wp(M) 



Ai/) 



w, 



p (6, y) Wp(6, j,) Wpl{yi b) + Wp{h) y) w ^ll$ 



<T a } 



Using the similar argument as above, we obtain (1191) . ■ 

Proof of Theorem 14.31 We let e = j* for a large integer N > 0. 
Using Lebesgue dominated convergence theorem, continuity and strong Markov 
property of X , we have 



E x {e-*»* ;X aa >x + a} = E^e" 9 *- ; r+ +Q < <7 Q } 

AT-1 



I + "-1,X + (i+l)E 






JV->oo 



i=0 



= exp lim 



AT-1 



N- 



E K+«{ 



_ gj4 *+»e-a,« + (i+l)e + 



' '.T+(i+l)e 



exp 



;=o 



£+a 



W qj -i(x,u—a) 
Wq(x.u-a) 



(ill 



s(u) — s(x) W q ,\ (x,u— a) 
s f (x) W q (x,u— a) 



Here we used Proposition ^, ll in the last equality. Moreover, notice that E x {e qA * 
r+} = E I {e-9 A r°' DO ; X T - < z}. It follows that, 

E x {e- qB °;X aa e{x,x + a)}= f E x {e~ qA ^ a '°° ; X T -_ 6 da} 

Jx 



x+a a 

dz 



Ex{- 



-qAl "• 



x+a 



s'(u) du 



x (, , s(u)-s(x) W q]1 {x,u-a) \ 2 W q (x, U - a) 

\ + s'(x) W q {x,u-a) ) 



where we used Proposition 14.11 in the last equality. 
Proof of Theorem 14.71 Let us denote by 



.g:=inf{t>0 : X t = X e J. 



Then we have 



E l = / hY t >y} at + / i {Yt>y} dt : = E q / + sy. 

J0 Jg 

Below we compute the Laplace transforms of E^ 1 and E q ' 2 conditioning on X Bq . 
More specifically, for m > x, we let e = ln j^ £ - for a large N > 0. Then we have 
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that 



E.je-P^' 1 ; X eq G dm} 



dm 



d_ 
"dh 



r + 



E x {exp I - p / l{Y t>y } di ) ; X 6g > to + /i} 
/i=o V -/o 

AT-l 



TT r / f T .+ C+i). \ , 

00 ;=n \ Jo / 



x — L-dP™(t + , , < e„) 



AT-1 



cxp( lim ^pE B+te {e" p/o "" t ' ' e ^<*+»-»> dt - ?T *% + i><} - 1] ) . 

V N—t-oc *• — ' J 



(j>q{m) 



exp 



m Wg, 2 (M-2/,u) + W g ,i(u,w- y) 



<t> q+P (it-v) 



- W qA (u -y,u) + W q (u, u~y) fe* ( ( ;_^ 



</>g"(m) 



■ dw 



The fourth equality follows from Corollary 14. 2 1 On the other hand, 
E x {e-^' 2 |X e =m} 



E m {exp -p 



— lim ■ 

e'->0+ 



*-{X t <m-y} ' 



dt h e «< T m+e'} 



= lim 



Pm(e q < T+ +£ ,) 



- m (e g < T m+e , ) 



^ T ™~V <e i AT rn+e>) E r p f«« l iXt<m _ y} dt < + 1 



1+ lim 



m (e g < T m+e , ) 

. Pm(e,<r+ +e ,) 



^^{[e-'/o- i<**<»-,>d*_i] . e? < T + } 



(37) 
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To get the limit in ([57]). we use Corollary 3.4 of [H? to proceed as 



lim 

:'-s-04 



= lim 



- rn{T m _ y < e q A T m+e , 



»( e 9 < T , 



■E m - v {[e- p ^ 1 < Jt «<»-»> dt -l]; 



e„ < t, 



m+e 



J 



my^q ^ 'm+e' 



W q {m + e',m) 
W q (m + e',m — y) 
0+(m) 
0q (m + e') 

W q (m + e', to — 3/) Y P ; — s'(m — y) + W 9) i(m — y,m + e') 



p + q 



W q {m + e',m- y) ^^-g + W, fl ( m - y, to + e') 



4>+(m-y) 
(j>q{m + e) 



(to) ^'( m ^y)+ i+p^M-^C" 1 -!/) 



•ftg + p (™-y) 



9# M W g ,i (m - y, m) + W g (m, ro - y) V p , ( " ^ 

It follows that 

E x {e-^' 2 |X e9 =m} 

Wg,2(m-y,m) + W 9 ,i(m,m - y) - ^s'(to) 



<# +p '("*-y) 



<+p(™-^) <#M 



W g>1 (m - y.ro) + W,(m,ro - y) *<£lZ-vl 
The proof is complete after integration with respect to to. i 



#j (m) 
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